The spatial range of the Kondo effect: a numerical analysis 
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The spatial length of the Kondo screening is still a controversial issue related to Kondo physics. While renor- 
malization group and Bethe Anzats solutions have provided detailed information about the thermodynamics of 
magnetic impurities, they are insufficient to study the effect on the surrounding electrons, i.e., the spatial range 
of the correlations created by the Kondo effect between the localized magnetic moment and the conduction elec- 
trons. The objective of this work is to present a quantitative way of measuring the extension of these correlations 
by studying their effect directly on the local density of states (LDOS) at arbitrary distances from the impurity. 
The numerical techniques used, the Embedded Cluster Approximation, the Finite U Slave Bosons, and Numeri- 
cal Renormalization Group, calculate the Green functions in real space. With this information, one can calculate 
how the local density of states away from the impurity is modified by its presence, below and above the Kondo 
temperature, and then estimate the range of the disturbances in the non-interacting Fermi sea due to the Kondo 
effect, and how it changes with the Kondo temperature Tk. The results obtained agree with results obtained 
through spin-spin correlations, showing that the LDOS captures the phenomenology of the Kondo cloud as well. 
To the best of our knowledge, it is the first time that the LDOS is used to estimate the extension of the Kondo 
cloud. 

PACS numbers: 73.23.Hk, 72.15.Qm, 73.63.Kv 



I. INTRODUCTION 



The physics of isolated impurities inside a Fermi sea has 
received considerable attention since it was experimentally 
shown that nano-systems composed by quantum dots pos- 
sess Kondo phenomena, very clearly reflected in its transport 
properties. 1 One signature of this effect is a narrow resonance, 
at the Fermi energy, in the local density of states (LDOS) 
of the impurity, with a width of the order of a characteristic 
energy, the so-called Kondo temperature, Tk- The transport 
properties of a nanoscopic structure in this regime are sub- 
stantially affected by the Kondo resonance, as it creates an 
extra channel at the Fermi level through which the electrons 
can propagate. The energy fcsTk is also associated with an 
antiferromagnetic correlation between the impurity and the 
conduction electron spins in its neighborhood, favoring the 
emergence of a singlet ground state. These spins, localized 
in the impurity's vicinity, constitute a screening cloud of the 
localized impurity spin, known as the Kondo cloud. While 
most of the physics involved in this important effect is by now 
well established, the nature, structure, and extension of the 
Kondo cloud, and even its existence, is still, to some extent, 
controversial^ 3 Theoretically, it is thought to be a crucial in- 
gredient in helping to understand, for instance, the interaction 
between two nearby impurities, when one of them is sitting 
within the region of influence of the Kondo cloud of the other. 
From the experimental point of view, although it is thought 
that the extension of this characteristic cloud can reach very 
large values? 2 - the properties of a system of impurities in met- 
als seem to depend linear on the impurity concentration. This 
seems to indicate that the impurities do not see each other, al- 
though, based on the expected Kondo cloud extension, they 
should. Moreover, there has not been any clear experimental 



evidence of it existence, with the exception of the electronic 
conductance measurements in quantum corrals,^ For instance, 
in an ellipsoidal quantum corral, a Kondo peak produced by 
a magnetic atom located at one focus of the ellipse has an 
experimentally detectable spectral response in the other fo- 
cus. This indicates a very peculiarly structured Kondo cloud, 
which, through the use of an Scanning Tunneling Microscope 
(STM), can be experimentally analyzed. 

The Kondo cloud length can be estimated by consider- 
ing that the mean-life of the Kondo quasi-particles are re- 
lated to the time-scale tr ~ ft/fcsTk. Assuming that these 
quasi-particles propagate with the Fermi velocity vp, then the 
Kondo screening length can be related to the quantity 2 



Rk 



Tivf 



(1) 



Obviously, since all electrons whose energies fall within the 
Kondo peak will participate in the formation of the Kondo 
cloud, the quantity vp is not well defined. Moreover, one 
may expect that the quasi-particles do not propagate with the 
bulk vp, but with a renormalized v*, given by the presence of 
the impurity. From Heavy Fermion theory, we can estimate 
v* = kp/m*, where m* is the effective mass of the quasi- 
particle? 4 . Therefore, Eq. [TJ should provide inaccurate results 
for the screening length Rk- However, it should be expected 
to give the correct dependence with Tk and some plausible 
order of magnitude for its length. 

From the theoretical point of view, this problem has been 
analyzed using different approaches i 6 - 7 - 8 ' 9 The study of spin 
properties, through the local susceptibility or the spatial spin 
correlation function, has given significant contributions to the 
understanding of this phenomenon. 610 More closely related 
to our approach, the analysis of the conductance of a quantum 
dot embedded in a finite wire* 2 - or the persistent currents in a 
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finite ringji using renormalization arguments or Density Ma- 
trix Renormalization Group calculations, respectively, were 
proposed as a way of determining the Kondo cloud, as well. 

More recently, a variational approach was proposed to 
study the propagation, from the impurity, of the local hole- 
density. 1 ' In this work, it was possible to show that, in two 
and three dimensions, the extension of the Kondo cloud is of 
the order of a few Fermi wave-lengths only, due to angular dis- 
persion effects, such that Rk does not play a significant role 
in the physics of a system of impurities in either of these di- 
mensions. This seems to explain the situation from the exper- 
imental point of view (as mentioned above), and thus the irrel- 
evance of the Kondo cloud in most of the real systems stud- 
ied. However, for one-dimensional systems^-!, the impurity- 
impurity interaction should be determined by the Kondo cloud 
length i?K- This will have important consequences to the con- 
ductance properties, and therefore will have implications to 
the design of quantum dot integrated nanoscopic systems. 

A. Holzner et ai, using DMRG, have calculated the spin- 
spin correlations involved in the formation of the Kondo cloud 
in a one-dimensional system and found the dependence of 
the range of the Kondo cloud with the Kondo temeperature 
to agree with Eq. QJ£ 

Motivated by this situation, we study the Kondo cloud in 
a one-dimensional system, focusing our attention on its elec- 
tronic properties. The study of the propagation of the Kondo 
resonance, located at the vicinity of the Fermi energy, will 
shine light, for instance, into the transport properties of a 
quantum dot connected to leads where the distance from the 
dot to an STM tip is changed in a controlled and continu- 
ous wayJ^ This will be experimentally similar to the trans- 
port properties studies of a system formed by a magnetic atom 
located in one focus of an elliptical quantum corral, as men- 
tioned above, and can be experimentally implemented for a 
quantum dot connected into an infinite wire. 

In the present work, we discuss the spatial behavior of 
the Kondo cloud by alternative means in an infinite one- 
dimensional system. Indeed, to estimate the cloud range, 
three different numerical techniques are used to track the ef- 
fects of the impurity over the LDOS far away from the im- 
purity. These effects are calculated above and below the 
Kondo temperature Tk, and their difference is used as a finger 
print of the extension of the Kondo cloud. The calculations 
are carried out using the Embedded Cluster Approximation 
(ECA), 14 the Finite-U Slave Bosons Mean-Field approxima- 
tion (FUSBMF), 16 and the Numerical Renormalization Group 
(NRG) method.ii 

The paper is organized as follows: In the next section 
(Sec. HI), we present the model used and the methods to solve 
it. In Sec. [HI] we briefly describe the behavior of the LDOS, in 
real space, within the metal lead and define the function used 
to estimate Rk- In Sec. IIVI the numerical results calculated 
using ECA, FUSBMF and NRG are discussed and compared. 
Finally, in the last Section, we present our conclusions. 
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FIG. 1: (Color online) Schematics of the system used to determine 
the spatial extension of the Kondo cloud. The exactly diagonalized 
(ED) cluster (with a variable number of sites L+l, where the QD is 
numbered as site zero) is indicated by the dashed box. Note that the 
(integer) index < N < oo runs over all sites inside and outside the 
cluster. By applying the numerical methods employed in this work, 
the dressed (interacting) Green function for all sites in the cluster is 
obtained (note that for FUSBMF and NRG, the cluster can be con- 
sidered as just the impurity - see end of Sec. Ill BB for details). A 
simple Equation of Motion procedure allows the calculation of the 
GF in any site of the semi-chain outside the cluster. With that, the 
effect of the Kondo screening effect over the non-interacting Fermi 
sea, the so-called Kondo cloud, can be probed at an arbitrary distance 
from the impurity. 



H. SYSTEM AND NUMERICAL METHODS 

In this work, we analyze a system composed by one Ander- 
son impurity (representing, for example, a quantum dot (QD) 
or an add-atom in a metal surface) coupled by a matrix ele- 
ment t' to a band (modeled by a semi-infinite non-interacting 
chain - from now on referred to as a semi-chain). This system 
is shown schematically in Fig.Q] This figure presents the sys- 
tem using the terminology appropriate for ECA, where a finite 
cluster has to be defined. However, as described in more de- 
tail below, we will show that, to calculate the Green functions 
(GF) outside this cluster, when doing FUSBMF or NRG, the 
same terminology can be used, although there is no equivalent 
cluster definition in FUSBMF or NRG. 

As shown in Fig.Q] N is an (integer) index that numbers the 
sites from zero to infinity, being the impurity (green circle), or 
QD, located at site N = 0. The letter L is not an index and its 
meaning, related to ECA, but extended to the other methods, 
is explained below in Sec III A fl 

The total Hamiltonian reads, 



H T - #imp + ^band + ^hybrid' 



(2) 



with 



ndcr^da (3) 

fT (7 
OO 

#band = 1 ( C Na C N+la + 4f +la C N<T ) (4) 

^hybrid = t'^iclvda +c\ a c 0a ), 



(5) 



where cl a creates an electron at the impurity, c ! Na . creates 
an electron at the site N of the lead, and n^a = cl a co a 
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is the number operator at the impurity. The first two terms 
of Ht represent the Hamiltonian of the impurity and the 
non-interacting band, respectively, and the last term is the 
hybridization between them. An important quantity for 
this system is the broadening of the impurity level V = 
27ri' 2 pi eac j(£'F), where Pi eal }(EF) is the LDOS of the first 
site in the semi-chain at the Fermi energy Ep. Note 

In this paper, we will concentrate our attention on the 
electron-hole symmetric point (V g = —U/2), although the 
results can be generalized to an arbitrary value of gate poten- 
tial. As mentioned in the introduction, we want to estimate the 
extension of the Kondo cloud, and its dependence on U/T, by 
analyzing the LDOS calculated through the local GF. To cal- 
culate the GF, we will use the ECA, the FUSBMF, and the 
NRG methods, which are briefly described next. 

At the end of this section we will describe how the LDOS 
in a site N far away from the impurity is calculated using 
the Equation of Motion (EOM) method. Note that the EOM 
method described in Sec. IIIBI does not depend on the method 
used to calculate G dd ; any of the three methods described be- 
low provide essentially the same kind of input for the EOM 
procedure. 

A. Numerical methods 

1. Embedded Cluster Approximation method 

The ECA method has been developed to treat localized im- 
purity systems consisting of a many-body interacting region 
weakly coupled to non-interacting conduction bands. The ap- 
proach is based on the idea that the many-body effects of the 
impurity are local in character (the Kondo cloud, for instance). 
With this in mind, we proceed in three steps: first, out of the 
complete system (the impurity plus a non-interacting band, 
described by a tight-binding Hamiltonian), one isolates a clus- 
ter consisting of the impurity plus their L nearest neighboring 
sites in the tight-binding semi-chain. This cluster, with a vari- 
able size L + 1, as it includes the impurity, is shown in Fig.Q] 
by dashed lines. The first site outside the cluster is labeled 
N = L + 1 (remember that N is an index, as opposed to a 
number of sites). Most of the many-body effects are expected 
to be confined to this cluster. 

The second step of the method consists in exactly diagonal- 
izing the cluster, using, for example, the Lanczos method^ 
and calculate all the GFs. Finally, in the third step, the clus- 
ter is embedded into the rest of the tight-binding semi-chain 
using a Dyson equation^ 

Being gij a cluster-GF that propagates an electron from 
site i to j and gL+i the GF of the first site out of the cluster 
(N = L + 1), the Dyson equation to calculate a dressed (by 
the presence of the semi-chain), GF for sites inside the cluster 
can be written as, 

Gis = gi,i + gi,L t Gl+i,{, (6) 
Gl+i,i = 9l+i t Gh,i ; (7) 

Note that the hopping parameter t in Eq. ^corresponds to the 
broken link shown by a dashed line in Fig. Q] This matrix 



element has the same value as all the other hopping parame- 
ters within the chain. Properties like conductance through the 
impurity and its LDOS, for example, can be obtained solving 
this set of equations (for more details, see Ref.flil). 

2. Finite- U slave bosons mean-field approximation 

The slave boson mean field is a method proposed origi- 
nally to treat the problem when the Coulomb repulsion U is 
the larger quantity. The double occupancy is excluded from 
the Hilbert space with the help of projectors-bosons opera- 
tors. After taking a mean field in the boson operators, the 
many-body Hamiltonian is mapped into an effective one-body 
Hamiltonian that can be solved exactly J£ 

The FUSBMF approach is an extension of the usual slave 
boson mean field in order to treat problems with finite 
The first step is to enlarge the Hilbert space, by introducing 
a set of slave boson operators e, p a and d, and replacing the 
creation and annihilation (dj CT ) operators in the Hamilto- 
nian by d\z\, and z a di a , respectively. Following Kotliar and 
Rukenstein^ the operator z takes the formic 

to = [l-d^d-pip^^p^+pU] 

x[l-e*e-plp„] 1 ' 2 . (8) 

Notice that the bosonic operators d and d) do not carry an 
spin index. The enlarged Hilbert space is then restricted to the 
physically meaningful subspace by imposing the constraints 

P = eJe + J2 vli>a + d f d -1 = (9) 

£7 

and 

Qa = ri da - pIp<t - d'd = 0. (10) 

Both constraints are included into the Hamiltonian through 
Lagrange multipliers A^ 1 ' and Xa ■ The constraint described 
by Eq. [9] force the dots to have empty, single, or double oc- 
cupancy only, while the constraint of Eq. [10] relates the bo- 
son with the fermion occupancies. In the mean-field approx- 
imation, we replace the boson operators e, p a and d (and 
their Hermitian conjugates) by their thermodynamical expec- 
tation values e = (e) = (et), p a = (p a ) = (p^) and 
d = (d) = (<r ). These expectation values, plus the Lagrange 
multipliers, constitute a set of parameters to be determined 
by minimizing the total energy (H). In principle, it is nec- 
essary a set of seven self-consistent parameters. Once again, 
as in the infinite U case, the problem was reduced to a one- 
body Hamiltonian whose energy can be minimized easily. The 
quantity we need to calculate is the Green function at the im- 
purity, around the Fermi level. Thus, 

G' M ={{z^M)), (ID 

which is the propagator that carries the correct weight of the 
Kondo resonance. 

In Sec. Ill Bl it will be shown how to calculate the GF in the 
lead's sites using G dd as an input. 
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3. Numerical renormalization group approach 

The NRG method was originally proposed by K. G. Wilson 
to study magnetic impurity problems^ Initially, it was ap- 
plied to the Kondo Hamiltonian, and later extended to the An- 
derson model^- 8 It can be shown that for these two models, at 
low temperatures, the states close to the Fermi level (i.e. with 
the lowest energy contribution) are the most relevant. There- 
fore, perturbation theories are not the most adequate approach 
to these problems. As a brief description of the method (a full 
detailed description can be found in Refs. Il7lll9l) . we present 
the two main steps in the implementation of the method.— 
The first one consists in sampling the energy interval of the 
conduction band by a set of logarithmically decreasing en- 
ergy intervals [xflf, x jv-i], defined by x^ = ±DA~ N , where 
A is the discretization parameter and D is the half-width of 
the conduction band. Then, from each interval, only one rep- 
resentative energy value is kept (chosen according to a well 
defined criterion, see Bulla in Ref. 12 for details). The total 
number of representative energies, one from each interval, re- 
sult in the set of discrete energies that couples to the impurity. 
After these two basic steps, the total Hamiltonian is mapped 
into a semi-infinite chain, commonly known as Wilson-chain, 
where each site of the chain corresponds to an energy scale in 
the logarithmically discretized conduction band, with the im- 
purity sitting at its first site. It is important to notice that the t n 
couplings, between adjacent sites n and n + 1, decrease, away 
from the impurity, as A~ Tl//2 . The final form for the Hamilto- 
nian in the NRG framework is 



H = lim A-^- 1 '/ 2 H 



N- 



N 



(12) 



where, 



H N = A^- 1 )/ 2 



coa + n.c. 



JV-1 



^ ] e n c l la -C na + ^ ] t n {c\ l<J C n J r i a + h.C. 



(13) 



where d a annihilates an electron with spin a at the impurity, 
and c na annihilates one at site n in the semi-infinite chain 
(indexed from n = to N). 

Note that an explicit analytical expression for t n in Eq. [13] 
cannot be obtained for a band of arbitrary shape. For the 
present problem, where a semi-elliptical band is used, we are 
forced to calculate the t n numerically. 20 The hoppings t„ that 
define the Wilson-chain must not be confused with the ma- 
trix elements t of the real space chain, shown in Fig Q] The 
elements t n correspond to the band obtained after the loga- 
rithmic discretization of the real space chain. It can be shown 
that when A — ► 1, the hoppings t n — ► tM 

The second important step consists in solving numerically 
the resulting Hamiltonian given by EqQ~2] To this end, we start 
with a system consisting of the isolated impurity, described by 
the Hamiltonian H[ mp . Then, the subsequent sites are added 
one by one. This procedure generates a sequence of Hamilto- 
nians Hn, which are solved as follows: at a given iteration N 



the Hamiltonian Hn is diagonalized numerically. The eigen- 
vectors and the corresponding eigenvalues are obtained. Next, 
a new site N + 1 is added. This is done by enlarging the cur- 
rent Hilbert space (associated to iteration N) through a ten- 
sorial product of its elements with the states of the site being 
added in the next iteration. This process results in an exponen- 
tial growth of the dimension of the Hilbert space of successive 
iterations. Due to computational constraints, it is necessary to 
truncate the Hilbert space at each iteration, after it reaches a 
certain size. The NRG truncation criterion is to keep only the 
A/ lowest energy states of Hm (typically, M = 1000), and 
neglect the higher energy spectrum. 

The process of adding a single site to is repeated un- 
til the system reaches the strong coupling fixed point. When 
this fixed point is reached, Hn and -ffjv+2 have the same 
eigenvalues^ 

The sequence of iterations described above can be thought 
of as a renormalization group (RG) process. Adding one site 
to the chain, and obtaining the new low-energy spectra, can be 
understood as an RG transformation 7Z that maps the Hamil- 
tonian Hn into a new Hamiltonian -ffjv+i = JZ(Hn), which 
has the same form as H^. Once the fixed points are obtained, 
the static and dynamic properties, as well as temperature ef- 
fects, can be calculatedJi^ In particular we are interested in 
the local GF at the impurity. 

At this point it is worth to remind the reader that the in- 
formation about the high energy dynamics is not accurately 
taken into account, since the high-energy spectra is partially 
neglected after the truncation. 

All the NRG data presented in this work was calculated 
with A = 2.5 and keeping the M = 1000 lower energy states 
in each iteration. To calculate the LDOS at the impurity, the 
delta functions were broadened using Logarithmic-Gaussians 
with a b = 0.6 factor (see Ref. HI). 

Finally, we want to stress, once again, the difference be- 
tween the real space semi-infinite chain and the Wilson-chain. 
The first one, shown in FigQ] has all the hopping terms equal 
to t. The Wilson-chain is just used to calculate the GF at 
the impurity, and it is obtained after the discretization of the 
real space chain. With the impurity propagator, obtained from 
NRG, the LDOS at any site of the real space chain can be 
calculated, as explained below. 



B. LDOS away from the impurity: Equation of Motion 

In this section, we will explain how to calculate the LDOS 
at any site of the semi-chain, which models the electron reser- 
voir. Note that, for all numerical methods used in this work, 
once the dressed GF is known at the impurity (and, in the 
case of ECA, for all the other sites of the cluster), a proce- 
dure based on the construction of a Dyson equation, through 
the use of a sequence of Equation of Motion, can yield the 
dressed GF for any site in the tight-binding semi-chain, no 
matter how far away from the impurity. This can be most eas- 
ily understood in the case of ECA, as this idea is built into the 
very core of the method. Indeed, ECA allows us to calculate 
not just the LDOS in all the sites of the cluster, but also in 
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all the sites in the rest of the semi-infinite tight-binding chain 
used to represent the lead. An important fact that we want 
to remark regarding ECA is that the embedding procedure re- 
sults in a feedback of the leads into the central region, but also 
reciprocally. The physics under study does not need to be re- 
stricted to entirely occurring within the exactly solved region. 
I.e., many-body effects taken into account exactly inside the 
ED cluster are propagated, by the Dyson equation, into the 
electron reservoir (the semi-chain), which now does not have 
anymore the LDOS of a non-interacting system. It is impor- 
tant to remark, as will be clearly explained shortly, that the 
change of the LDOS in the semi-chain from tight-binding to 
many-body lies at the core of the method used in this work 
to estimate the range of the Kondo cloud. One added benefit 
of the procedure to be described below is that the physics of 
the Kondo effect at the impurity (the Kondo resonance) does 
not need to be calculated with ECA for the EOM procedure to 
work. In the present work, it is also calculated with FUSBMF 
and NRG. 

To calculate the dressed propagators at any site of the semi- 
chain, we write down the EOM of the local propagators at a 
site M, In the case of ECA, the site M must be outside the 
cluster, i.e., M > L + 1 (see Fig. [T). This restriction does 
not apply to FUSBMF or NRG, where the equivalent to the 
ECA cluster can be considered to be just the impurity. A brief 
description of the EOM method can be found in Ref. To 
simplify the notation, in what follows we will ignore the spin 
index a. The set of equations to solve, in order to calculate 
G m ,m, is given by 



G 



Gm,m 
Gm+i,m 
Gm-i,m 

Af+l.Af-l 



,9o + got Gm-i,m + 9ot Gm+i,m, (14) 
g sc t Gm,m, (15) 

Gj\/,M-1 

got Gm-x,m-x + got Gm+i,ai-i, (16) 
g sc tGM.M-i, (17) 



where go = 1/uj is the atomic GF at site M and g sc is the 
bare propagator for the rest of the semi-chain starting at the 
site M + 1 and is given by, 



± Vu 



At 2 



2t 2 



(18) 



In Eq[l6l we used explicitly the equivalence between 
Gm-\,m and Gm,m-i- This is only valid if the hopping pa- 
rameters t are real (e.g., no magnetic field inside the chain, al- 
though a more general EOM, involving a magnetic field, can 
also be found. 

Solving this set of equations, we obtain, 



go 



G 



M,M 



got got 
1 - got 2 g s 



Gm- 



l.M-l 



1 - g t 2 g s 



(19) 



Note that the GF at the site M can be calculated after the GF 
at M — 1. Note that the equation above clearly indicates that, 
to calculate the dressed GF at site M, the only many-body 
information needed is the dressed GF at site M — 1. This fact 
automatically defines a procedure to find the propagator at any 



site in the semi-chain. Note that Eq. [19] is defined for a site, 
within the semi-chain (i.e. M > 2), that is connected to both 
adjacent sites by a matrix element t. Thus, this still leaves 
us with the task of calculating G\ t \. To calculate the correct 
propagator at site 1, we have to rewrite Eqs.[l4]to[T7] in order 
to obtain Gi t i as a function of 6*0,0, without overlooking that 
the hopping between sites and 1 is t', not t. 

For FUSBMF and NRG, we start with the GF calculated at 
the impurity, i.e., Go,o- Using the EOM method, we calculate 
the propagator at the first site of the chain, G\.\. Then, using 
Eq. [19] the propagator Gm,m can be calculated at any site. 
The procedure for the ECA method is slightly different, as 
in ECA all the dressed propagators inside the ED cluster are 
calculated already within the method. Therefore, in ECA, the 
EOM procedure starts at site L + 1 (see Fig. [O, using Eq.[T9l 
where Gl.l is an input from the ECA calculations. In that 
case, there is no special procedure to calculate G^i. 

Once Gm,m is calculated for the desired site M, the LDOS 
can be calculated as, 



Qm{u) = — Im[G M ,Af] 

7T 



(20) 



As we are using the same procedure to find the LDOS 
away from the impurity for three very diverse numerical meth- 
ods (ECA, FUSBMF, and NRG), some explanation about 
the adopted terminology is necessary, so that the same term, 
with slightly different meanings, can be unambiguously used 
throughout the manuscript. As explained in Fig. Q] in ECA, 
cluster means a variable size finite group of sites (including 
the impurity), which is exactly diagonalized and embedded 
(as explained above). In this manuscript, the ECA cluster 
contains up to L + 1 = 10 sites (i.e., the impurity plus up 
to L = 9 tight-binding sites). For sites N > L + 1, the LDOS 
will be found through the EOM method, as described above. 
On the other hand, an FUSBMF or an NRG cluster, given the 
very nature of both methods, contains just the impurity itself 
(therefore, L = 0, see Fig.[T|). Because of that, there is a slight 
difference to the application of the EOM method to these last 
two methods, viz., Gi,i has to be calculated first, and then all 
the other Gi.i are calculated by using Eq. [I9]in sequence, as 
explained above. 



III. LOCAL DENSITY OF STATES WITHIN THE METAL 
HOST 

Using the FUSBMF approximation, we can obtain analyt- 
ically self- consistent expressions for the GF, and then the 
LDOS. The local GF is also obtainable within NRG, i.e., the 
LDOS at the impurity can be found with very good accu- 
racy, and, as explained in Sec III Bl the GF (and therefore the 
LDOS) can be calculated within the semi-chain that models 
the non-interacting band. The same is valid for ECA, despite 
the distinctions drawn above between ECA, in the one hand, 
and FUSBMF and NRG, in the other hand. 

In Fig. 12 we show the LDOS for several different sites 
within the non-interacting semi-chain (N > 0), calculated by 
ECA. The parameters used are U = t and T = Q.lt, and 
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FIG. 3: (Color online) (a) LDOS for site N = 0, where the impurity 
is located (see Fig. [TJ, calculated by ECA. Solid line shows LDOS 
for V g = —(7/2, U — t, and t' = 0.3t, while in dashed lines, re- 
sults are shown for T > Tk and the same parameter values for V a , 
U and t' . An imaginary part v — 0.001 was used to regularize the 
LDOS. The dashed (black) curve displays an example of the normal- 
ized Lorentzian used to convolute the data in Eq. [22] (b) LDOS for 
the first and third sites of a non-interacting semi-chain with band- 
width D — At. (c) Same as (b), but now for the second and fourth 
sites. 



FIG. 2: (Color online) (a-f) The LDOS for sites N = 1,2, 3, 4, 50, 
and 100 (see Fig. [TJ, calculated by ECA for a cluster with L = 3 
(note that the impurity is located at TV = 0). Solid (black) curves 
show the LDOS at T = for V g = -(7/2, U = t, and t' = 0.3t, 
while dashed (red) curves show results for T > Tk and the same 
values for V g , U, and t'. An imaginary part v — 0.001 was used to 
regularize the LDOS. Note that the LDOS's for the first three panels 
were calculated using ED, in contrast to the ones for the last three 
panels, which were calculated using the EOM method described in 
the text. The LDOS for the first three sites (N = 1, 2, and 3) was 
also calculated with EOM and, as expected, there was a very good 
agreement between the results obtained with the two different meth- 
ods. 



the gate potential is set at the electron-hole symmetric point 
(V g = —(7/2). The LDOS for sites outside the cluster were 
obtained through the procedure described in Sec lIIBI The 
LDOS at zero-temperature is shown in solid lines. In dashed 
lines, for comparison, we show the LDOS when the impurity 
is out of the Kondo regime (T > Tk). In order to obtain a 
solution for T > Tk, we use the Hubbard-I approximation, 23 
which artificially eliminates the spin correlations between the 
impurity and the leads. This approximation is equivalent to 
performing an ECA calculation where the cluster contains just 
the impurity (L = in Fig. Q] i.e., the atomic solution). To 
understand the results in Fig. [2 it is instructive to analyze the 
results shown in Fig. [3j where, in panel (a) it is shown the 
LDOS at the impurity, for T = [(red) solid line) and for 
T > Tk [(blue) dotted curve], and in panels (b) and (c) the 



LDOS of the first four edge sites of an isolated non-interacting 
semi-chain. The results shown in Fig. [2] [panels (a) to (d)] dis- 
play essentially the hybridization between the LDOS of the 
impurity (panel (a) in Fig. and the LDOS of the sites in 
the semi-chain (panels (b) and (c) in Fig. [3j, once the im- 
purity is coupled to the semi-chain. The solid (black) curve 
shows the LDOS for T = 0, while the dashed (red) curve 
shows the LDOS for T > Tk, at the first four sites in the 
semi-chain, after it couples to the impurity. This hybridiza- 
tion can be described in a simple way: a peak in the LDOS 
of the impurity [Fig. EJa)], centered at co p , will generate ei- 
ther a resonance or an anti-resonance (at u) p ) in the LDOS of 
a semi-chain site when the impurity couples (hybridizes) to 
the semi-chain. On the one hand, a resonance (a peak) will re- 
sult if the semi-chain's LDOS, in one specific site, vanishes at 
u) p . On the other hand, an anti-resonance (a dip) results when 
the site's LDOS at cu p is finite. This resonance/anti-resonance 
site to site oscillation effect in the LDOS will have important 
consequences in the next section. Our interest is to be able 
to distinguish the effect caused over the semi-chain's LDOS, 
far away from the impurity, by the presence (T < Tk) or 
absence (T > Tk) of a resonance at the Fermi energy (the 
Kondo peak) in the LDOS at the impurity [compare the solid 
and dotted curves in Fig.[3ja)]. The extent to which this hy- 
bridization effect can spread away from the impurity will be 
used as a measure of the extent of the Kondo cloud. 

Figure [2] clearly displays this kind of hybridization effect, 
as described above. Indeed, if one concentrates the attention 
on the features close to the Fermi energy (ui = 0) in the differ- 
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FIG. 4: (Color online) (a) Main panel: Absolute value of the cloud 
extension function, |T(N)|, calculated by ECA as a function of the 
semi-chain site N for several different cluster sizes (L = 1 to 9, in 
steps of 2 - see legend). The parameter values are U = t and Y = 
Q.lt. Note that L is the number of nearest neighboring sites to the 
impurity within the cluster solved exactly in ECA. The curves shown 
are the data obtained through Eq.|22] Fits to these data (not shown), 
for each value of L, using Eq.[23] result in exactly the same curves 
as the ones shown. It is clear that the data for each different cluster 
size decays exponentially with N (notice the logarithmic scale in the 
vertical axis). The open (red) triangles curve shows the extrapolation 
of |.F(N)| to the thermodynamical limit, as describe in detail in the 
text. The inset shows the site to site oscillations of -F(N), as well as 
the exponential decay, now with a linear scale for the vertical axis, 
(b) ECA extrapolation to the thermodynamical limit (dashed (red) 
curve), as described in the text, of the Rk data (solid (black) dots) 
obtained from the fittings for each different L curve shown in panel 
(a), as explained in the text. 



(a) ECA (U=t) 




FIG. 5: (Color online) Cloud extension function F(N) as a function 
of the site in the semi-chain, for several values of T. (a) Extrapola- 
tion to the thermodynamical limit for each V value and (b) FUSBMF 
results. Note that, as V increases, the reach of the distortion in the 
LDOS of a given site N (-Rk), produced by the impurity, is reduced, 
reflecting the shorter range of the Kondo effect. Note also that the 
dependence of |_F(A)| with N for all T values shown is perfectly 
linear. 



ture the change from the solid to the dashed curve, say, in site 
50 [panel (e)], as that occurring in the LDOS away from the 
impurity when the temperature is lowered below Ik. Panel 
(f), where there is very little difference between both curves, 
shows that the impurity, in a Kondo regime, has a spatially 
limited influence over the Fermi sea. It is one of the aims 
of this paper to understand how this influence depends on the 
sole energy scale of the Kondo effect, i.e., the Kondo temper- 
ature Tk- 



ent panels in Fig. [21 one sees that the difference between the 
LDOS curves below Tk (solid) and above Tk (dashed) is quite 
marked, and owes its origin to the presence of the Kondo peak 
at the impurity below Tk- By using the Lorentzian shown in 
Fig. [3] (dashed line) to restrict ones attention to the immediate 
neighborhood of the Fermi energy, by convoluting it with the 
difference between the solid and dashed curves in Fig.|2j one 
expects to extract the essence of the influence of the impurity, 
when in the Kondo regime, over the Fermi sea. One can pic- 



IV. NUMERICAL RESULTS 

As mentioned in Sec. QJ considering that the width of the 
Kondo resonance in the LDOS at the impurity, A, is propor- 
tional to Tk, we expect that, 
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FIG. 6: (Color online) Range (Rk) of the Kondo cloud as a func- 
tion of A. The open dots show the results obtained for Rk through 
the fitting of curves like the ones in Fig. [5] for various values of T 
(therefore, different A, i.e., (the width of the resonance at the Fermi 
level A), (a) ECA, (b) FUSBMF, and (c) NRG results. A fixed pa- 
rameter U = t was used for the three methods, and the parameter T 
was changed in order to obtain different A values. The solid (red) 
line in each panel shows the interpolation of an 1/A function, as 
expected for Rk vs Tk (which is proportional to A). The value of 
the proportionality factor for each method is shown in the respective 
label. 
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FIG. 7: (Color online) ECA results for the size of the Kondo cloud 
vs. 1/A for U = OM, U = t and U = 1.5t. Note that, for the 
interval of values of U and T shown, Rk is exactly a linear function 
of 1/A, and does not depend on the values of U and V used. 



In this section, we will estimate the screening length Rk by 
evaluating the distortion in the LDOS, produced by the Kondo 
resonance at the impurity, in sites N arbitrarily far away from 
the impurity. The distortion produce in site N will be quan- 
tified by the absolute value of the function F(N), as defined 

by 



/oo 
-oo 



(22) 



where g§ (g^ K ) is the local density of states at the site N in 
the Kondo regime (out of Kondo), and La is the Lorentzian 
distribution of width A [dashed curve in Fig. [3ja)]. Note that 
the distance to the impurity is given by r — Na, where a is 
the lattice parameter. It is important at this point to remind the 
reader that all the calculations here were done at V g = —U/2 
(i.e., at the particle-hole symmetric point). To evaluate g^ K , 
we applied the EOM procedure to the Hubbard-I solution for 
the impurity's GF (where the dotted line in Fig.[3ta) shows the 
negative of its imaginary part). 

Note that, for a site far away from the impurity, where the 
many-body effects are not important, g§ must be equal to 
g^ K , thus F(N) sa 0. This can be seen in the last panel of 
Fig. |2] F(N) will be used to find a length scale beyond which 
the presence of the impurity is not relevant any more. We will 
call _F(N) the cloud extension function. 

As already mentioned, to calculate g^, we will use ECA, 
FUSBMF, and NRG. 

In Fig. Ufa), it is shown the absolute value of the cloud ex- 
tension function, |F(N)|, calculated with ECA, as a function 
of N, for U = t,T = O.lt, and for several values of L (note 
that the vertical axis has a logarithmic scale). The scatter plots 
show the data obtained from Eq.|22]for some selected values 
of N, for different cluster sizes L (see legend). We find that 
the behavior of |F(N)|, for all values of L used, is a decay- 
ing exponential (the correlation factor, when fitting the curves 
with an exponential, was exactly 1 for all values of L): 



\F(N)\=A cM-N/Rk) 



(23) 



where Rk marks the distance from the impurity where the 
value of |-F(N)| has fallen by 1/e, in comparison to its value 
at the impurity (N = 0). The solid lines in panel (a) show 
the result of fitting each set of data with Eq. [23] In the inset 
of Fig. |4ja), it can be observed that F(N) oscillates between 
positive and negative values for successive NM- This is a direct 
consequence of the resonance/anti -resonance oscillation in the 
LDOS discussed in the previous section. The extrapolation of 
Rk to an 'infinite' ECA cluster (1/L — > 0) is explained in 
detail next. 

As mentioned already, each curve in Fig.HJa) (from L = 1 
to L = 9) is fitted using Eq.|23] Therefore, a value for Rk(L) 
is found for each cluster size used in ECA (1 < L < 9). 
FigureHfb) shows the values obtained this way for Rk(L) as 
a function of 1 /L (solid (black) dots), for the different val- 
ues of L used in panel (a). A fitting of these results by a 
quadratic polynomial is also shown (dashed (red) curve). The 
intercept of the dashed (red) curve with the vertical axis pro- 
vides an extrapolation of Rk to the thermodynamical limit 
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Rk(L — > oo). This extrapolated value of Rk can then be used 
in Eq.[23]to obtain the thermodynamical limit for the |F(7V)| 
curve, which, for Y/t = 0.1, is the open (red) triangles in 
Fig.ffia)^ 

Obviously, Eq. [23] has two free parameters, viz., Aq and 
Rk- Although the vertical axis in Fig. |4fa) is logarithmic, 
making it difficult to judge the convergence of Aq, it is true 
that A converges with L faster than Rk (note that, in accor- 
dance with Eq.[23] A is the y-intercept and Rk is the negative 
of the inverse of the slope of the curves for different cluster 
sizes). The values of Aq for different values of L where ob- 
tained from the fitting, of each data set in Fig.lUa), done with 
Eq. [23] The dependence of Aq on the model parameters will 
be discussed below^ 

At this point, it is important to note that, as the cluster used 
in the FUSBMF and NRG calculations has a fixed size (L = 
0), there is no extrapolation to be done to find |.F(iV)| for both 
methods. There is only a fitting to Eq.[23]to find Rk and Aq. 

Now that we have clarified how the thermodynamical limit 
value for Rk is found for each method, we want to show how 
it varies with the model parameters. Figure[5la) shows the ex- 
trapolated (ECA) \F(N)\ curves, for U = 1.0 and different 
values of Y (from 0.0625 (solid (black) curve) to 0.4 (dashed 
(magenta) curve). From the data, it is clear that Rk (the nega- 
tive of the inverse of the slope) decreases with increasing val- 
ues of r. In panel (b), the corresponding |.F(iV)| curves ob- 
tained with FUSBMF are shown, for comparison. The over- 
all agreement between both methods is quite good. As ex- 
pected, we observe that the size of the Kondo cloud (measured 
through the cloud extension function) increases with U/Y. We 
can understand this behavior by noting that, as Y increases 
(with a fixed U), 2k also increases, and Rk, as predicted by 
Eq.Q] decays. 

Figure [6] shows the results obtained [open dots in panels 

(a) , (b) and (c)] by extrapolating Rk from |F(7Y)| for differ- 
ent values of Y (at V g = -U/2), using ECA (a), FUSBMF 

(b) , and NRG (c). These results are plotted as a function of A 
(which is taken as the full-width at half-height of the Kondo 
peak for each different value of Y. We observe that the depen- 
dence of the Kondo length Rk with A satisfies the relation- 
ship given by Eq. [21] (as A, the width of the Kondo peak, is 
proportional to Xk). To emphasize that, each set of data (ob- 
tained by the three different methods) was fitted by a function 
oc 1/A (see the solid (red) line in each panel). It is important 
to stress that the proportionality coefficient between Rk and 
1/A obtained by all the three different methods is very sim- 
ilar, i.e., Rk ~ 2.0/A. The proportionallity factor in Eq.l2TI 
are 1.874 for ECA, 1.964 for FUSBMF, and 2.102 for NRG. 
While this factor is similar for ECA and FUSBMF, there is 
a 10% difference between ECA and NRG. We believe that 
this difference comes from the parameter b used in NRG to 
broaden the logarithmic-Gaussian functions in the LDOS, as 
the value of A obtained by NRG is very sensitive to the choice 
of this arbitrary parameter. 

Figure [7] shows ECA results of Rk v.s. 1/A for three dif- 
ferent values of U (0.5, 1.0, and 1.5). As shown in the fig- 
ure, the functional form Rk oc 1/ A is valid for the intervals 
of U and A (and therefore Y) inside which the calculations 
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FIG. 8: (Color online) Parameter Aq as a function of A (in log- 
scale) calculated by ECA and FUSBMF. Different U where used for 
ECA. Note that, when A is small, all curves coincide, as expected, 
in view of the universal behavior characteristic of a Kondo system. 
Increasing A, we enter in the mixed valence regime, and curves for 
different values of U start to diverge. 



were done. More importantly, as all the curves collapse to a 
single line, the proportionality coefficient is also independent 
of these intervals. This indicate that the 'propagation' of the 
Kondo effect into the leads (which is essentially measured by 
Rk) depends only on the weight of the Kondo resonance at 
the Fermi level (measured by A). 

Figure [8] shows, for ECA and FUSBMF, the parameter Ao 
as a function of A. We can see that, the ECA and FUSBMF 
curves agree quite well, for the same value U = 1.0. Ad- 
ditional calculations, with different U values (U = 0.5 and 
U = 1.5), were done just with ECA. For these additional re- 
sults, one sees that the curves start to differ from each other 
for large A, but agree for small values (A < 0.04). The agree- 
ment for small A can be easily understood if one takes in ac- 
count the universal behavior of the Kondo effect, in the sense 
that it is determined by a single energy scale, the Kondo tem- 
perature Tk- Therefore, Aq (as Rk) does not depend on ei- 
ther U or r independently, but on their ratio (U/Y), at least 
until the system enters the mixed valence regime, at higher A 
(equivalent to Y). Notice that the curve for lower U (U = 0.5) 
starts to diverge from the other two at a lower value of A 
(proportional to Y), while the opposite occurs for the larger-C/ 

27 

curve.-' 



V. CONCLUSIONS 

Using, for the first time, the LDOS within the metal lead, 
we have estimated the effective length, in real space, of the 
effect of the many-body correlations originating at the impu- 
rity site. For the Kondo effect, we defined a cloud extension 
function F(N) in order to estimate when an electron located 
at a site N, away from the impurity, is not affected anymore 
by its presence. 
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Indeed, we have used the electronic properties reflected in 
the LDOS function (charge spectra) of the one-dimensional 
metallic lead, to study the spatial propagation of the Kondo 
effect away from the magnetic impurity. The length of the 
Kondo cloud, Rk, has been defined in the literature to be the 
extension of the spin-screening-cloud, formed by the conduc- 
tion electrons, in the vicinity of the impurity. From this point 
of view, it is essentially the spatial size of a magnetic property, 
as it is associated to the spin-spin correlations between the 
local impurity and the conduction electron spins However, 
here we claim that, as far as the Kondo cloud is concerned, the 
charge spectra counterpart of the Kondo physics (defined as 
the effect of the impurity's Kondo peak over the LDOS of the 
leads) is equivalent to its magnetic expression, as they are both 
manifestations of the same physical phenomenon. Moreover, 
we claim that \F(N)\, being dependent just on the LDOS, 
is easier to calculate and measure than spin-correlation-based 
functions. 

As far as transport properties are concerned, if the Kondo 
effect is thought to be the way in which two or more QDs can 
interfere, the relevant way of studying the Kondo cloud is by 
analyzing the effect of the impurity over the LDOS of the rest 
of the system. 

In order to study the Kondo effect spatial propagation, we 
define what we call a cloud extension function, denoted F(N) 
(see Eq. |22| |. It measures, in an interval of width Tk around 
the Fermi energy, the distortion of the LDOS, at site N, cre- 
ated by a Kondo impurity sitting at the origin. We evaluated 
this function using three totally different formalisms, ECA, 
FUSBMF, and NRG, and obtained almost identical results for 
the variation of Rk with the parameters of the system. The 
fact that three different formalisms provide the same physi- 



cal description makes this study quite robust and reliable. We 
demonstrate, as well, that the length of the Kondo cloud is 
controlled by the unique, scaling invariant, relevant parameter 
of the Kondo effect, the Kondo temperature Tk- These results 
permit a very accurate determination of the functional form 
of Rk(Tk), in agreement with intuitive ideas, summarized in 
Eq.ID 

Finally, it is important to emphasize that the measurement 
of spin-spin correlations between different sites, in a real STM 
experiment, is difficult to perform, as the use of two different 
STM tips, at the same time, is required.— On the other hand, 
the mapping of the Kondo cloud through the difference in the 
conductance, measured by an STM tip, at different points in 
a ID system, looks more feasible, as it has already been per- 
formed in metallic surfaces^ 
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from the semi-chain (by making t = 0). Our calculations were 
done both ways, i. e., Hubbard-I and t' — 0, and we found agree- 
ment between the two results for the extension of the Kondo cloud 
Rk- 

Note that, in ECA, A was calculated as the thermodynamical limit 
-1/L — > 0- of the width of the Kondo resonance at the impurity's 
LDOS. 

A second extrapolation method can also be used (although not 
as accurate as the one described in Sec lIVb : one can extrapolate 
the value of | ^(TV) | to the thermodynamical limit (1/L — > 0) 
for each value of N, and then obtain Rk by fitting the extrapo- 
lated |F(iV)| to Eq.f23] Although there is good qualitative agree- 
ment on the values obtained for Rk between the two extrapolation 
methods, we favor the first one, as it involves much less computa- 



tional effort and is therefore more accurate. 
The limit in which |F"(iV)| — > determines the order of magni- 
tude of the Kondo cloud. However, the actual values and func- 
tional variation of both Aq and Rk in Eq. [23] clearly indicate 
that Rk has a leading contribution in determining the size of the 
Kondo cloud, given the small variation of Ao in the region of pa- 
rameter space associated to the Kondo effect (compare the scales 
of the vertical axis in Figs. 6 and 8). 

Note that Ao was not calculated as |F(0)], but through the fittings 
discussed in Figs. 4 and 5. Therefore, its value, and its dependence 
with A, as seen in Fig. 8, do not represent the value or the varia- 
tion of the area of the Kondo peak. 
M. Ternes, private communication. 



